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Abstract 

The Kowalevski top in two constant fields is known as the unique profound example of an integrable Hamil- 
tonian system with three degrees of freedom not reducible to a family of systems in fewer dimensions. As the 
first approach to topological analysis of this system we find the critical set of the integral map; this set con- 
sists of the trajectories with number of frequencies less than three. We obtain the equations of the bifurcation 
diagram in R 3 . A correspondence to the Appelrot classes in the classical Kowalevski problem is established. 
The admissible regions for the values of the first integrals are found in the form of some inequalities of general 
character and boundary conditions for the induced diagrams on energy levels. 

1 Introduction 

During the last 20 years the integrable case of S. Kowalevski [1] has received several generalizations. Among them 
a special place is given to the case [5] of rotation about a fixed point of a heavy electrically charged gyrostat in 
gravitational and electric fields of force. For a rigid body without gyrostatic effects the corresponding equations 
were first considered in [3] and interpreted as the equations of motion of a massive magnet subject to the gravity 
force and constant magnetic force fields. The mathematical model of superposition of such fields is referred to as 
^ ' two constant fields [4] . 

The case [2] does not have any explicit groups of symmetry and therefore provides an illustration of a physically 
realizable system with three degrees of freedom not admitting any obvious reduction to a family of systems with 
two degrees of freedom. The phase topology of irreducible systems has not been studied yet. The theory of n- 
dimensional integrable systems originated in [5] has not been further developed due to the absence, at that moment, 
of non-trivial natural examples. 

The result [2] succeeded some previous publications dealing with rigid bodies and gyrostats satisfying the 
conditions of Kowalevski type: I.V. Komarov [5] has proved the complete integrability of the Kowalevski gyrostat 
in gravity force field by finding the first generalization of the Kowalevski integral K; the corresponding integral 
for the rigid body in two constant fields was pointed out by O.I. Bogoyavlensky [3j; this integral was upgraded to 
the case of gyrostat by H. Yehia [7J. Yet the analog of the Kowalevski case for two constant fields had not been 
considered integrable until A.G. Reyman and M.A. Semcnov-Tian-Shansky [2] found the Lax representation with 
spectral parameter; this immediately led to the new integral generalizing the square of momentum integral for 
axially symmetric force fields. 

Later, in a joint publication with A.I. Bobenko [4], the authors of [2] presented algebraic foundations for 
the integrability of multidimensional Kowalevski gyrostats and described a viable way of explicit integration using 
finite-band technique. For two constant fields this integration was never fulfilled. 

This paper starts the investigation of three-dimensional phase topology of a rigid body of Kowalevski type in 
two constant fields. 
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2 Preliminaries 

Consider a rigid body with fixed point O. Choose a trihedral at O rotating along with the body and refer to it 
all vector and tensor objects. Denote by eie2e3 the canonical unit basis in R 3 ; then the moving trihedral itself is 
represented as Oeie 2 e 3 . 

Constant field is a force field inducing the rotating moment about O of the form 

r x a (2.1) 
*REGULAR AND CHAOTIC DYNAMICS, V. 10, N 4, 2005, pp. 381-398, DOI: 10.1070/RD2005v010n04ABEH000321 



1 



with constant vector r and with a corresponding to some physical vector fixed in inertial space; r points from O 
to the center of application of the field, a is the field's intensity. 

For two constant fields the rotating moment is ri x a + r 2 x (3. ft can be represented as (|2 . 1 [) if either 
ri x r 2 = or a x [3 = 0. In the sequel we suppose that 

ri x r 2 ^ 0, axj3/0. (2.2) 

Two constant fields satisfying (|2.2p are said to be independent. 
Introduce some notation. 

Let L(n, k) be the space of n x fc-matrices. Put L(k) = L(k, k). 

Identify R 6 = R 3 x R 3 with 1,(3, 2) by the isomorphism j that joins two columns 

A = j(ai,a 2 ) = ||ai a 2 || € £(3,2), ai,a 2 € R 3 . 

For the inverse map, we write 

j-\A) = ( Cl (A),c 2 (A)) e R 3 x R 3 , A e L(3,2). 

If A, B G L(3, 2), a S R 3 , by definition, put 

2 

/lxB = ^c,(4)xc,(B)eR 3 ; a x A = j(a x Ci(A), a x c 2 (A)) e L(3, 2). (2.3) 
i=i 

Lemma 1. Let A e 50(3), _D e GL(2, R), a e R 3 ; A, 5 e L(3, 2). Then 

A(A xB) = (AA) x (AB); (AD- 1 ) x (BD T ) = AxB: 
A(a x A) = (Aa) x (AA); a x (AD) = (a x A)£>. 

The proof is by direct calculation. 

Denote by I the inertia tensor of the body at O and by us the angular velocity. Using the notation (|2.3p we 
write the Euler - Poisson equations of motion in the form 

la;" = Ia> X a> + A X Z7, IT =-u>xU. (2.4) 

Here A = j(ri, r 2 ) is a constant matrix, ?7 = j(cx,(3). The phase space of (|2.4j) is {(w, C/)} = R 3 x L(3, 2). 
In fact, [/ in (|2.4p is restricted by geometric integrals; that is, for some constant symmetric C € L(2) 

U T U = C. (2.5) 

Let O represent the set (|2 . 5[) in i(3, 2). In order to emphasize the C-dependence, we write O = O(C). 
Let S = (I, A, C). Denote by X s the vector field on R 3 x 0(C) corresponding to the system (|2.4p . 
Associate to A e 50(3), D e GL(2,R) the linear automorphisms *(A, D) and %p(A,D) of R 3 x L(3,2) and 
L(3) x L(3,2) x L(2) 

y(A,D)(u>,U) = (Au,,AUD T ), 

ip(A, D)(I, A, C) = (AIA T ,AAD-\DCD T ). [ > 

It is easy to see that ([22]) and (jUSJ) imply *(A,D)(R 3 x O(O)) = R 3 x 0(DCD T ). Using Lemma ED we 
obtain the following statement. 

Lemma 2. For each (A,D) E 50(3) x GL(2,R), we have 

9{A,D).(X s (v)) = X mtD){s) (V(A,D)(v)), veTL 3 x 0(C). 

Thus, any two problems of rigid body dynamics in two constant fields (for short, RBD-problems) determined 
by the sets of parameters 5 and i/)(A, D)(S) are completely equivalent. 

Let us call an RBD-problem canonical if the centers of application of forces lie on the first two axes of the 
moving trihedral at unit distance from O and the intensities of the forces are orthogonal to each other. 

Proposition 1. For each RBD-problem with independent forces there exists an equivalent canonical problem. 
Moreover, in both equivalent problems the centers of application of forces belong to the same plane in the body 
containing the fixed point. 
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Proof. Let the RBD-problem determined by the set of parameters S = (I, A, C) satisfy (|2.2|) . This means that the 
symmetric matrices A* = (A T A) -1 and C are positively definite. 

According to the well-known fact from linear algebra, A* and C can be reduced, respectively, to the identity 
matrix and to a diagonal matrix via the same conjugation operator 

DA*D T = E, DCD T = diag{a 2 , b 2 }, D G GL(2, R), a, b G R+. 

Then Ci(AD~ 1 ) and c 2 (AD _1 ) form an orthonormal pair in R 3 . There exists A G 50(3) such that 
Aci(AD _1 ) = ei(i = 1,2). The first statement is obtained by applying Lemma [2] with the previously chosen 
A, D to the initial vector field Xs- 

To finish the proof, notice that the transformation A h-+ AD^ 1 preserves the plane spanning a (A), 02(A). 
The matrix A in Q2.6|) represents the change of the moving trihedral. Therefore, if a G R 3 represents some physical 
vector in the initial problem, then Aa is the same vector with respect to the body in the equivalent problem. □ 

Remark. The fact that any RBD-problem can be reduced to the problem with one of the pairs 1*1, r2 or 
ol,(3 orthonormal is known from [8]. Simultaneous orthogonalization of both pairs crucially simplifies calculations 
below. 

It follows from Proposition [1] that, without loss of generality, for independent forces we may suppose 

ri = e l5 r 2 = e 2 , (2.7) 
a-a = a 2 , /3-/3 = 6 2 , at- = 0. (2.8) 

Change, if necessary, the order of ei, e 2 (with simultaneous change of the direction of 63) to obtain a b > 0. 

Consider a dynamically symmetric top in two constant fields with the centers of application of forces in the 
equatorial plane of its inertia ellipsoid. Choose a moving trihedral such that Oe^ is the symmetry axis. Then the 
inertia tensor I becomes diagonal. Let a — b. For any 9 G SO(2) denote by £ 5*0(3) the corresponding rotation 
of R 3 about Oe 3 . Take in $nj§ A = 9, D = 6. Under the conditions (|2"7T|) . (|2~%]l tp = Id and * becomes the 
symmetry group. The system (|2.4p has the cyclic integral Iuj • (a 2 e3 — a x (3) pointed out in [7] for the analog 
of the Kowalevski case. Therefore it is possible to reduce such an RBD-problem to a family of systems with two 
degrees of freedom. 

Let us call an RBD-problem irreducible if for its canonical representation (|2.7[) . (12. 8[) the following inequality 

holds 

a > b > 0. (2.9) 

The following statements are needed in the future; they also reveal some features of a wide class of RBD- 
problems. 

Lemma 3. In an irreducible RBD-problem, the body has exactly four equilibria. 

Proof. The set of singular points of (|2.4|) is defined by uj — 0, A x U = 0. For the equivalent canonical problem 
with (f27T|) 

e! x a + e 2 x/3 = 0. (2.10) 

Then the four vectors in (|2.10| are parallel to the same plane and |ei x ct| = |e 2 x (3\. With (|2.8|) . (|2.9p this 
equality yields 

a = ±ae 1 , [3 = ±6e 2 . (2.11) 

From mechanical point of view, the result is absolutely clear: none of the orthogonal forces with unequal 
intensities and "orthonormal" centers of application can produce a non-zero moment at an equilibrium. □ 

Lemma 4. Let an irreducible RBD-problem in its canonical form have the diagonal inertia tensor I. Then the 
body has the following families of periodic motions of pendulum type 



(2.12) 
(2.13) 



a — ±aei, uj = ip'ei, (3 = fe(e 2 cost/? — e3 sin ip), 
2ip" = —bsimp; 

(3 = ±fee 2 , uj = ip'e 2 , a = a(ei cos ip -\- e 3 simp) , 
2ip" = —a shop; 

a x (3 = ±a6e 3 , uj = <p'e3, 
a = a(ei cosip — e 2 sin 1^2), f3 = ±6(ei simp + e 2 cos^), (2-14) 
<p>" = — (a ± b) sin ip. 

The proof is obvious. Note that in the case considered the pointed out families are the only motions with 
constant direction of the angular velocity. In particular, the body in two independent constant fields does not have 
any uniform rotations. 
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3 Critical set of the Kowalevski top in two constant fields 



Suppose that the irreducible RBD-problem has a diagonal inertia tensor with principal moments of inertia satisfying 
the ratio 2:2:1; then we obtain the integrable case of the Kowalevski top in two constant fields. By an appropriate 
choice of measurement units, we present equations (I2.4[) in scalar form 



2uj\ = uu 2 u 3 + fa, = -uixu>3 - a 3 , u> 3 = a 2 - fa, 
a\ = a 2 uj 3 - ct 3 u> 2 , fa = fau 3 - fau 2 , 
a 2 = a 3 ui - aiuj 3 , fa 2 = fawi - fau> 3 , 
a 3 = aiuj 2 — a 2 uji, fa 3 = faoj 2 — fauj\. 



(3.1) 



The phase space is P 6 = R 3 x O, where O c R 3 x R 3 is defined by O is diffeomorphic to SO(3). 

The complete set of first integrals in involution on P e consists of the energy integral H, the generalized 
Kowalevski integral K [3] , and the integral G found in [5] : 

H = uj\+u)1 + ^1- (ai+/3 2 ), 

K = {w{ -ujj+ai- fa) 2 + {2u x u 2 +a 2 + Pi) 2 , 

1 1 (3 2) 

G = -(2aiwi + 2a 2 u 2 + a 3 oj 3 ) 2 + -(2faw x + 2fauj 2 + facu 3 ) 2 + v ' ' 

+ ^uj 3 (2-f 1 uj 1 + 2-f 2 uj 2 + j 3 lj 3 ) - b 2 ai - a 2 fa. 

Here we denote by 7, the components of 7 = a x f3 relative to the moving basis. 
Introduce the integral map 

J = G x K x H : P 6 ^ R 3 . (3.3) 

Let a C P 6 be the set of critical points of J. By definition, the bifurcation diagram of J is the subset £ C R 3 
over which J fails to be locally trivial; £ determines the cases when the topological type of the integral manifolds 

J c = J-\c), c= (g,k,h) eR 3 (3.4) 

changes. Finding the critical set a and the bifurcation diagram is the necessary step in the topological analysis of 
the problem as a whole. 

It follows from the Liouville - Arnold theorem that for c ^ S the manifold (|3.4p . if not empty, is a union of 
three-dimensional tori. The considered Hamiltonian system is non-degenerate (at least for sufficiently small values 
of b); then the trajectories on such a torus are quasi-periodic with three almost everywhere independent frequencies. 
The critical set a is invariant under the phase flow and consists of trajectories with number of frequencies less than 
three. These trajectories are called critical motions. For a generic value c £ S the set J c D a consists of two- 
dimensional tori. The dynamical system induced on the union of such tori for c in some open subset in S is a 
Hamiltonian system with two degrees of freedom. Vice versa, let M be a submanifold of P 6 , dimAf = 4, and 
suppose that the induced system on M is Hamiltonian. Then, obviously, M C o . This speculation gives a useful 
tool to find out whether a common level of functions consists of critical points of J. 

Lemma 5. Consider a system of equations 

h = 0, f 2 = (3.5) 

on a domain W open in P 6 . Let X be the vector field on P e corresponding to \3. 1}) and M C W defined by iS.5\) . 
Suppose 

(i) fi an d f'2 o,re smooth functions independent on M ; 
(n) Xfi = 0, Xf 2 = on M; 

(Hi) the Poisson bracket {/1, f 2 } is non-zero almost everywhere on M. 
Then M consists of critical points of the map J. 

Proof. Conditions (i), (ii) imply that M is a smooth four-dimensional manifold invariant under the restriction of 
the phase flow to the open set W . Condition (iii) means that the closed 2-form induced on M by the symplectic 
structure on P 6 is almost everywhere non-degenerate. Thus the flow on M is almost everywhere Hamiltonian with 
two degrees of freedom. It inherits the property of complete integrability. Then almost all its integral manifolds 
consist of two-dimensional tori and necessarily lie in a. Since M is closed in W and a is closed in P 6 , we conclude 
that M Ca. □ 

Two systems of the type p. 51) are known. The first one was pointed out in [3]. It is the zero level of the 
integral K. The condition K = leads to two independent equations defining the smooth four-dimensional manifold 
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97t C a. It is shown in [9] that the 2-form induced on 971 by the symplectic structure on P 6 is degenerate on the 
surface of codimension 1. 

The second critical subset 01 C a was found in [10] in the form of a system of two equations satisfying the 
conditions of Lemma [5j The functions in these equations have essential singularities at the points 

ai = 2 , a 2 = -Pi- (3.6) 

The set 01 was investigated in [TT] . It was shown that 01 is the set of critical points of some smooth function F 
on P 6 . Then 01 is stratified by the rank of Hesse's matrix of F and fails to be a smooth four-dimensional manifold at 
some points of the set (|3.6p . In particular, it cannot be defined by any global system of two independent equations. 
In this case the induced 2-form also has degenerate points even in the smooth part of 01. 

The following result completes the description of the critical set a by adding a new invariant subset O C P 6 ; D 
is almost everywhere a smooth four-dimensional manifold. Note that the sets OJt, 01 and D have pairwise nonempty 
intersections corresponding to bifurcations of critical integral manifolds of the induced " almost Hamiltonian" sys- 
tems with two degrees of freedom. 

Let us introduce the following notation 

p 2 = a 2 + b 2 , r 2 = a 2 — b 2 ; 
6 = a-L - 2 , 6 = «2 + 01, Vi = «i + P2, i] 2 =a 2 - Pi- 

Theorem 1. The set of critical points of the integral map &3.3\) consists of the following subsets in P 6 : 

1 ) the set 9Jt defined by the system 

Zi = 0, Z 2 = (3.7) 

with 

Zi =w?-c^ + &, Za = 2wiwa + £a; (3.8) 

2) the set 01 defined by the system 

Fi = 0, F 2 = 0, £ + i\ + (3.9) 



with 



and by the system 



Fi = (£i 2 + £f)w 3 - 2[(fiwi + t 2 io 2 )a 3 + - ^ 2 )0 3 ], 

F 2 = " £ 2 2 )(2^ 2 + 6) ~ 2&£ 2 (w? - lu 2 + Ci), 



(3.10) 



6=6 = 0, a 3 = ±r, 3 = 0, V 2 + r, 2 = 2(p 2 - r 2 ), 
(uj 2 + uj 2 )(a 3 Lj 3 + rjiLOi + r\ 2 Lo 2 ) + r 2 uji = 0; 

3) the set O defined by the system 

p x = 0, R 2 = (3.12) 

with 

Ri = (a 3 u 2 - 3 loi)lj 3 + 2£iUiuj 2 - &(lj 2 - J 2 ) + r) 2 {u>\ + wf), 

R 2 = {a 3 uoi + 3 uj 2 )lu 2 + [a 2 + 2 + ^(lu 2 - lu 2 ) + 2f 2 WiWa+ (3.13) 
+ Vi( u i + u 2 )]lo 3 + 2[6 (a 3 wi - /3 3 w 2 ) + 6 ("3^2 + A^i)]- 

Proof. Introduce the change of variables [10 (i 2 = —1) 



!El=6+«6 5 £2=6-^2, 
2/1 = ??i + "/2, 2/2 = »7i - w?2, 
zi = a 3 + i/3 3 , z 2 = a 3 -i/?3, 
ioi = wi + iu>2, W2 = — iu> 2 , w 3 = lo 3 . 



(3.14) 



The system (|3.1j) takes the form 



x'j = -I1W3 + ZlWl, x' 2 = X 2 W 3 - Z 2 W 2 , 

y'l = -yiw 3 + z 2 wi, y' 2 = y 2 w 3 - ziw 2 , 
2z[ = xiw 2 - y 2 wi, 2z' 2 = -x 2 wi+yiw 2 , 
2w[ = -(wiw 3 + zi), 2w' 2 =w 2 w 3 + z 2 , 2w' 3 = y 2 -yi. 



(3.15) 



Here the prime stands for d/d(it). 
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Denote by V 9 the subspace of C 9 denned by (|3.14ll . On V 9 , equations (|2.8p of the phase space P 6 become 



z\ + xiy 2 = r 2 , z 2 + x 2 yi=r 2 



x 1 x 2 + y 1 y 2 + 2z 1 z 2 = 2p 2 . (3.16) 
By virtue of (|3 . 14[) and (|3 . 16[) the integrals (|3.2|) take the form 

H = + w i w 2 - ^(i/i + 3/2), 

^ !/ 2 x 2 1/ x (3-17) 

G = -(p - xxa; 2 )u; 3 + -(a^i^i + x 1 z 2 w 2 )w 3 + 

+ ^(i2««i + Uiw 2 )(y 2 wi + X1W2) - -^P 2 ivi + V2) + -^r 2 (xi + x 2 ). 

Let / be an arbitrary function on V 9 . For brevity, the term "critical point of /" will always mean a critical 
point of the restriction of / to P 6 . Similarly, df means the restriction of the differential of / to the set of vectors 
tangent to P 6 . 

While calculating critical points of various functions (in the above sense), it is convenient to avoid introducing 
Lagrange multipliers for the restrictions (|3.16[) . Notice that the following vector fields 



d 






X 2 — , A3 — - 

dw 2 c 


d 1 d 


Z2 ^~ 
ax 2 


+ Zl- TT^l^— 

2 dz\ 


d 




dxi 


+ Z2 Ti o^ 2 ^ - 

dm 2 dzi 


d 


d d 


Xl-x — 

ox\ 


dx 2 dyi 



d_ 

>w 3 



1 d 

2 V1 ^ 
1 d 

2 X2 dz 2 ' 
d 

dy 2 

are tangent to P 6 C V 9 and linearly independent at any point of P e . Then the set of critical points of / is defined 
by the system of equations 

Xtf = 0, X 2 f = 0, X 3 f = 0, (3.18) 
Yxf = 0, Y 2 f = 0, Y 3 f = Q. (3.19) 

1. Apply p.ISp and (|3.19p to / = K. Then a critical point of K satisfies either 

w\ + x x =0, w\ + x 2 = (3.20) 

or 

wi = w 2 — 0, zi — z 2 = 0. (3-21) 

The system (I3.20P coincides with (|3.7I) and the only invariant set generated by (|3.21[) consists of all points of 
the trajectories (|2.14p . Such points satisfy (|3.12p . 

2. Consider the regular points of K at which H and K are dependent. Applying (|3.18p to f = H + sK with 
Lagrange multiplier s we immediately obtain w 3 — 0. Then from (|3.15p we come to solutions (|2.12p . (|2.13p . Along 
the corresponding trajectories both conditions (|3.9p . p,12p are valid. 

3. We now assume that H and K are independent. Introduce the function with Lagrange multipliers r, s 

L = 2G + (t - p 2 )H + sK. 

The multiplier of G is non-zero by assumption. The term with p 2 is added for convenience. 
The set 00 C a of the points satisfying for some r, s the condition 

2dG + (r - p 2 )dH + sdK = (3.22) 

is preserved by the phase flow of (|3.15p . Applying the corresponding Lie derivative to (I3.22p gives 

r'dH + s'dK = 0. 

Since dH and dK are supposed to be linearly independent, on <jq we obtain 

t = 0, s = 0. (3.23) 
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Hence r, s are partial integrals of motion on the invariant surface do- 
Equations (|3.18|) with / = L give 



x 2 ziw 3 + x 2 y 2 wi + (r - z 1 z 2 )w2 + 2stUi(u| + x 2 ) = 0, , ^ 

xiz 2 w 3 + (r - ziz 2 )wi + xiy\w 2 + 2sw 2 (w\ + x x ) = 0; 

(r — X1X2V3 + a^ziu"! + a;iz 2 w 2 = 0. (3.25) 

First consider the case (J376J . From (|3 . 1 6[) we come to the following values 

x x =x 2 = 0, z\ = z\ = r 2 , y lV2 = 2(p 2 - r 2 ). (3.26) 

Equations (|3.24p and (|3.25p hold if wiw 2 = or w 3 — 0. If either of these equalities takes place on some 
interval of time (and hence identically), then we obtain one of the solutions (|2.12p - (|2.14|) . 
Let wiw 2 7^ 0, w 3 7^ at some point satisfying (|3.26|) . Then (|3.24|) and (|3.25[) yield 

t = 0, s = r 2 /(2 Wl w 2 ). (3.27) 

Since z\ and z 2 are complex conjugates of each other, it follows from (|3.26| that they are real and equal. 
Denote their value by z = ±r. 

With (|3.26[) and (|3.2T|) the system (|3 . 1 9[) reduces to a single equation 

wiw 2 [2zw 3 + {w 2 y 1 + wiy 2 )] + r 2 (w 1 + w 2 ) = 0, (3.28) 

which corresponds to (I3.11|) . 



Note that (|3 . 28[) is obtained from (|3 . 9[) as p — y £ 2 + ^| tends to zero only after dividing by the maximal 
available power of p. Thus, at the points (13. 6|) the system (|3.9|) . without the assumption that p ^ 0, has extra 
solutions not belonging to a. 

Suppose x\x 2 7^ 0. The determinant of (|3.24p with respect to r, s equals 6 = 2{x\w 2 — x 2 w\). Let S = 
on some interval of time; calculating the derivatives of this identity in virtue of (|3.15p . we obtain one of the cases 
(|3~20|) . (|3~2Tj) . Therefore we may assume that 5^0. Then (|3~24|) implies 

1 

■[(x 2 z 1 wi - X!Z 2 w 2 )w 3 + x 2 y 2 w\ - xiyiw 2 ], (3.29) 



2{x\w\ — x 2 w\ 

t = z x z 2 H 5—^ 7T {[x 1 x 2 (z 2 w 1 - ziw 2 ) - wiw 2 (x 2 z 1 w 1 - xiz 2 w 2 )]w 3 - (3.30) 

X\W 2 — x 2 w{ 

- wiw 2 {x 2 y 2 w\ - xiyiw 2 ) + x^w^^yx - y 2 )}. 
Eliminating r from (|3 . 25[) and (|3 . 30[) we obtain Si = 0, where 

Si = [xiX 2 (z 2 Wl — ZlW 2 ) — WlW 2 (x 2 ZlWl — Xiz 2 w 2 )]w1 + 

+ [{xix 2 — ziz 2 )(x 2 wl — xiw 2 ) ~ wiw 2 (x 2 y 2 w\ — xiyiw 2 ) + 
+ xix 2 wiw 2 {yi - y 2 )]w 3 - (x 2 w\ - xiwl){x 2 ziwi + xiz 2 w 2 ). 

Next we solve (|3.25p for r and calculate the derivative r' in virtue of l|3.15p . According to (|3.23p we must 
have 5 2 = 0, where 

S 2 = (x 2 ziwi - xiz 2 w 2 )w 2 + (x 2 y 2 w\ - xiyiw 2 + x 2 z\ - xiz 2 )w 3 - 

- (yi - y2){x 2 ziwi + xiz 2 w 2 ). 

Notice that 

Si + wiw 2 S 2 = FiR. 

Here 

Fi = xix 2 w 3 - {x 2 ziwi + xiz 2 w 2 ) (3.31) 
is the first function from (|3.10p . The function 

R = {z 2 wi - ziw 2 )w 3 + x 2 w\ - xiw\ + wiw 2 {yi - y 2 ) (3.32) 

is a multiple of the first function from (|3.13p . precisely, R — 2iR\. Thus on the trajectories consisting of critical 
points, we have either F\ = or Ri = 0. Calculating the derivatives of these identities in virtue of (|3. 15[) we obtain 
(|3.9p and (|3.12p . respectively. Hence (|3.9p and (|3.12p provide necessary conditions for a point to belong to <jq. 

To prove sufficiency, it is enough to check (|3.19p . We avoid this technically complicated procedure and only 
notice that the systems (|3.9p and (|3.12p satisfy the assumptions of Lemma [5j □ 
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The phase topology of the induced system on OJl was studied in [9] . The system of invariant relations ()3.9() 
corresponds to that found in [10]. In the paper [11] the equations of motion on 0T are separated and the initial 
phase variables are expressed via two auxiliary variables, the latter being elliptic functions of time. The motions 
on 971 generalize those of the 1st Appelrot class (Delone class) of the Kowalevski problem [TJ]. As b tends to zero 
the motions on 01, as shown in [10] , convert to the so-called especially marvelous motions of the 2nd and 3rd classes 
of Appelrot [T2]. The set defined by the system (|3.12p was not pointed out earlier. 

To find the classical analog of the set D, put /3 = in (|3.13p . Then £x = r]i = a%, £2 = V2 = a 2 an d we obtain 

Rx = 2£uj 2 , R2 = 2£(uiu 3 + a 3 ), 

where 2£ — 2a,\U)\ + 2a2^2 + 0x30)3 is a constant of the momentum integral existing in the case of one force field. 
Therefore equations (|3.12p yield either £ — or 

u>2 = 0, + «3 = 0. (3.33) 

The condition £ = follows from the fact that when (3 = the integral G takes the value £ 2 . Equations ()3.33|) 
define especially marvelous motions of the 4th class of Appelrot [T2] . 



4 Bifurcation diagram 

Since all common levels of the first integrals (|3.2p are compact, the bifurcation diagram £ coincides with the set 
of critical values of the map (|3.3|) , that is, £ = J{u). 
Let 7 = \a X /3|. According to l|2.8[) . 7 = ab. 

Denote by A the region of existence of motions, that is, the set of c = (<?, k, h) 6 R 3 for which the integral 
manifolds (|3.4[) are not empty. 

Theorem 2. The bifurcation diagram of the map G x K x H is the intersection of A with the union of the surfaces 

T x : fc = 0; (4.1) 
T 2 : p 2 h- 2g + r 2 Vk = 0; (4.2) 
T 3 : p 2 h~ 2g-r 2 Vk = 0; (4.3) 

s e R\{0} (4.4) 



-Y 2 

k = 'is 2 -Ahs+p 2 + h 2 -^- 
2 3 

g = -s 3 + hs 2 + — 

s 



and the line segment 



T 5 : g = 7/1, k = p l - 27, h z < 4 7 . (4.5) 
In the parametric representation of the surface Y± the parameter s stands for a multiple root of the polynomial 

$( s ) = S A - 2hs 3 + [h 2 +p 2 - k)s 2 - 2gs + 7 2 . (4.6) 



Proof. 1. The equation of the surface (|4.1[) follows immediately from (|3.7[) , (|3.8p . and the expression of K in (|3 . 2|) . 
2. Relations (|4.2[) . (|4.3[) are equivalent to 

(p 2 h-2g) 2 -r 4 fc = 0. (4.7) 

Introduce the function 

F = (p 2 H - 2G) 2 - r 4 K. 

For X1X2 7^ denote 



tfi = t/ — K + zi), C/ 2 = J-K 2 +x 2 ) (C/ 2 = C/i). (4.8) 

" Si V ^2 



From representations (|3.17p and (|3.3ip we obtain 



p 2 H- 2G + r 2 VK = Ff + 2r 2 {lmU 1 ) 2 , 

p 2 H -2G -r 2 \fK = F? -2r 2 (KeU 1 ) 2 . 

2X\X2 
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Here v K is the principal square root of K. 

The equation of the zero level of F splits into two distinct equations 



F? + 4r 2 x 1 x 2 {lmU 1 ) 2 = 0, (4.10) 
Fl - 4r 2 xia;2(ReC/i) 2 = 0. (4.11) 



From (f3~10|) and ((48)) we have 



F 2 = ^^-(U? - U 2 ) = 2x x x 2 Im Ui Re Ul 
2i 

Thus the solutions of (|3.9p satisfy either l|4.10p or (|4.11[) and therefore lie on the zero level of the function F. 
The corresponding values of the first integrals satisfy (|4.7[) . 

From (|3.17p it follows that (|4.7[) holds for all points of the phase space such that X\X 2 =0 (regardless of their 
critical or regular nature). Hence (|4.7[) holds for the points (|3.11j) . 

3. Consider the system (|3.12[) . In terms of the variables (|3.14l) it is equivalent to the following equations: 

R = 0, R* = 0. (4.12) 

Here R is defined by (|3.32p and 

R* = (z 2 wi + ZiW 2 )w 2 + \x 2 w\ + x\w\ + w\w 2 {yi + y 2 ) + 2ziz 2 ]w3+ 
+ 2{x 2 ziWi + xiz 2 w 2 ). 

Notice that, after several differentiations in virtue of the system (|3. 15|) . the possibility z 2 wf — z\w\ = leads 
to the conditions (|3.2ip . that is, to the critical motions (|2.14p . Assuming (|3.21|) we obtain from (|3. 1 6[) 

x\x 2 = p 2 - 2q, yi2/2 = P 2 + 2q, (xi + x 2 )yxy 2 = r 2 (j/i + y 2 ) 

(q = ±7)- 

The corresponding values of the integrals (|3.17p are 

h = ^w 2 3 - + 2/ 2 ) > -Vp 2 + 2q, k=p 2 - 2q, g = qh. (4.14) 



If q = —7, then all of these values satisfy (J4T4J) with 

1, 



s = -\h± y/h 2 +4 7 ]. (4.15) 
Let q = 7. Then the values (|4.14p satisfy with 

s= hh± V/i 2 -47], (4.16) 

that is, only for the energy range h 2 47. For q = 7 and h 2 ^ 47 the values (|4.14p fill the segment (|4.5p . 

Consider the trajectories for which the equalities (|3.2ip do not hold identically. Express W3 from the first 
equation ijiTTZ)) : 

w 3 = \x 2 w\ - xiw 2 + w\w 2 (y\ - y 2 )]. (4-17) 

z 2 w\ - ziw 2 

Replacing u> 3 in (z 2 wi — ziw 2 ) 2 R* by (|4.17p . we obtain the expression 2wiW 2 Q (the resultant of (|3.32p . (|4.13p 
as polynomials in w 3 ), where Q is a non- homogeneous polynomial of third degree in wi and w 2 whose coefficients 
are polynomials in Xj, j/j and ^ of degree not greater than four. Since (|3.2ip is already excluded, the system p,12p 
is replaced by (|4.17p and the equation 

Q = 0. (4.18) 
We claim that in virtue of (|4.17p and (|4.18p . the values (|3.29p and (|3.30p satisfy the identities 

T-p 2 - 2s(s-H) = 0, 

(r - p 2 ) 2 + 4(p 2 - K)s 2 - 8Gs + {p A - r A ) = 0, (4.19) 
(t - p 2 )(2s - H) + 2(p 2 - K)s - 2G = 0. 

Here the calculation sequence is as follows. 
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We substitute (|3TT|) . (|3?29| . (f3T30|) . (f4TT7j) in the left-hand side of each equation ([419]) and multiply the result 
by the denominator, which is already supposed to be non-zero. The expression thereby obtained appears to be the 
product of some polynomial in variables (|3.14|) and the polynomial Q, which equals zero due to (|4.18|) . 

Replace in (|4.19[) the functions G,K,H by their constant values g,k,h and exclude r with the help of the 
first relation. The remaining two reduce to the form 



$(.s) = 0, d$(s)/ds = 0, 
where $ is the polynomial (|4.6p . Equations (|4.4|) are equivalent to (|4.20[) . 



(4.20) 
□ 



Remark. It is easy to see now that the relations (|4.ip - (|4.3[) turn into corresponding relations of the 1st, 
2nd, and 3rd classes of Appelrot as (3 tends to zero. Simultaneously, the polynomial (|4.6[) turns to stp(s), where 
tfi(s) is the Euler resolvent of the second polynomial of Kowalevski. This provides an alternative insight into the 
connection of the set T4 with the 4th Appelrot class of motions. The part of the segment defined by the 
inequality h 2 < 4j for the classical case (7 = 0) disappears. 



5 The region of existence of motions 

The results of the previous section are not complete until we find some conditions that give a criterion to establish 
whether a point of E = Ti U T2 U 1^ U T4 U Ts belongs to the region of existence of motions A = J(P 6 ) C R 3 . 
Three inequalities of general character can be obtained immediately from (13. 2| and (|4.9| : 



k > 0; 
ft ^ -(a + b); 



P h > 2g ■ 



r 2 ^fk. 



(5.1) 
(5.2) 

(5.3) 



the inequality obtained from (|5.3[) was established 



In case of the Kowalevski top in the gravity field (p 2 = 
by Appelrot [H] . 

To get more precise estimations for (g, k, ft) S A, restrict the problem to iso-energetic surfaces Eh 
P 6 : H(v) = ft}. Denote 

J h = GxK\ Eh :E h ^H 2 



{v e 



and let S/,, Eft, A^ be the cross-sections of E, E, A by the plane parallel to and height ft above the (g, fc)-plane. For 
any ft the set E/t is a bifurcation diagram of the map Jh- 

Notice that all sets Eh are compact. As proved in [13], they are connected as well. Therefore the values of 
any continuous function on Eh fill a bounded and connected segment. 

Let 

fe*(ft) =minif \E h , k*(h) = maxK \E h ; 
g*(h) = minG \E h , g*{h) = maxG \E h . 



(5.4) 



Then the rectangle 



n(ft) = {{k,g) : K(h) < A; < fe*(ft), g*(h) ^g^ g*(h)} 



cuts Eft out of E/j and we hope that this operation is not ambiguous. 

The following statements allow us to find explicitly the values (|5.4[) and give some more information about 
the sets r.nA. 

We are going to investigate various maps constructed of combinations of the first integrals G, H and, 
possibly, restricted to invariant submanifolds in P 6 . For each map /:(£—> R fe of this type we call a point c G R fe 
on admissible value if / _1 (c) 7^ 0. 

In Sections [3l H] we often referred to the motions (|2.12p - (|2.14ll . They will be also important in the sequel. 
Calculating the related values of G, K, H we obtain the sets in (ft, g)-plane and \ii in (ft, fc)-plane (i = 1, 6): 



Ai 


g = a 2 h + a(a 2 — b 2 ), 


Mi 


k = 


(ft + 2a) 2 , 
(ft - 2a) 2 , 
(ft + 2&) 2 , 
(ft-2&) 2 , 
(a-6) 2 , 
{a + b) 2 , 


ft ^ 


-{a- 


f b); 


A 2 


g = a 2 h — a(a 2 — 6 2 ), 




k = 


ft ^ 


a — I 




A 3 


g = b 2 h-b{a 2 -b 2 ), 




k = 


ft ^ 


-(0- 


f by, 


A 4 


g = b 2 h + b(a 2 -6 2 ), 


IH 


k = 


ft ^ 


-a-f 


b; 


A 5 


g = abh, 


M5 


k = 


ft ^ 


-(0- 


f 6); 


A 6 


g = —abh, 


Me 


k = 


ft ^ 


-a-f 


■6. 



The existing pairwise intersections of the first four sets in either group correspond to the equilibria (|2.11[) . 



10 



Recall that 97t = {K = 0} c P 6 and J(Wl) = Ti n A. Denote 

M = p 2 iJ - 2G : P 6 -> R 
and let = AfW = M\ m , The following result belongs to D.B. Zotev [9]. 

Proposition 2. (%) XTie function has three critical values hi — —2b, h-2 = 26 and h$ = 2a. In particular, 



mini? = -2b. 
wi 



(ii) The bifurcation diagram of — iJW x : 971 — > R 2 consists of the half-line 

m = 0, /i > -26 

and £/ie sei o/ solutions of the equation 



27m 4 + 4/i(/i 2 - 18p 2 )m 3 - 2[4p 2 /i 4 - (16p 4 + 15r 4 )/i 2 + 2p 2 (8p 4 - 9r 4 )]m 2 + 
+Ar i h[h 4 - Ap 2 h 2 + 2(2p 4 - 3r 4 )]m - r 8 [{h 2 - 2p 2 ) 2 - 4r 4 ] = 

in the quadrant {m ^ 0, h —26}. 

(m) XTie sei of admissible values of J' 1 ' is 

^ m ^ m (/i), ^ ^ -26, 

where mo(/i) stands for the greatest positive root of i5.6]) . which is considered as an equation in m. 

The bifurcation diagram of JW is shown in Fig. 1. The admissible values fill the shaded region. 



(5.5) 
(5.6) 





m 






m (h) 




-2b 


2b 


2« h 



Figure 1: The bifurcation diagram of i/W x MW 

The proof given in [9] is based on an ingenious change of variables on 971. Let us point out the relation between 
this result and Theorem [2] 
Let 

m=p 2 h-2g. (5.7) 

It follows from (|5.3p that m ^ on 971. This inequality explains (|5.5|) . Moreover, the line m = in the plane 
fc = is the intersection T\ n (T2 U F3) (in fact, along this line Ti and r 2 U are tangent to each other). 
The intersection Ti n F4 is defined by the system obtained from (|4.1[) . (|4.4p 

3s 4 - 4hs 3 + (p 2 + /i 2 )s 2 - 7 2 = 0, 



0. 



(5.8) 



s 4 — hs fi + gs — 7 

By virtue of the notation (|5 . T[) the left-hand side of (|5.6p becomes the resultant of the polynomials in 
with respect to s. Thus the set (|5.6[) corresponds to ri n Y4. 

Recall that 91 C P 6 is defined by 1(33)1 . (j3TT|l and J(9t) = (r 2 U T 3 ) n A. Let i/( 2 ) = G (2) = G\ m 

Introduce the map 

j( 2 ) =ij( 2 ) xG' 2 ' : 9T^R 2 . 



The following statement is proved in [14] 
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Proposition 3. (i) The bifurcation diagram of consists of the half-lines Ai, A 2 , A3, A4, the half-line 

1 



g = -p 2 h, h > -2b, 



and the curve 



2p 2 (p 2 h - 2g) 2 - 2r 4 h(p 2 h - 2g) + r 8 = 0, p 2 h ^ 2g. 
(ii) The admissible values of fill the region defined by the system of inequalities 

b 2 h - b{a 2 -b 2 )^g^ a 2 h + a(a 2 - b 2 ), h^-(a + b), 
2p 2 (p 2 h - 2g) 2 - 2r i h(p 2 h - 2g) + r 8 ^ 0. 

The bifurcation diagram of is shown in Fig. 2. The admissible values fill the shaded region. 




Figure 2: The bifurcation diagram of x G (2) 

Propositions [2] [3] completely define those parts of Ti, T2 and which correspond to real critical motions, 
that is, the sets ^nA and (r 2 U T 3 ) n A. 
Consider the map 

J (3) = H xK : P 6 -> R 2 . 

The critical set of is already found (see steps 1, 2 in the proof of Theorem [TJ. In addition to the manifold 
Wl it contains all pendulum motions (|2. 12|) - (|2.14|) . 

Proposition 4. (i) The bifurcation diagram of the map 

J(3) 

consists of the parabolic curves /13, [14, the 

half-lines n§, (Iq, and the half-line 

k = 0, h > -2b. (5.9) 



(ii) Let 



{h + 2b) 2 , -{a + b) < h < -2b 
0, h > -2b 

k*{h) = (h + 2a) 2 . 



k*(h) 



The admissible values of fill the region 

K(h) < fc < k*(h), h^ -(a + b). 



(5.10) 



(5.11) 



The inequality in (|5.9p follows from Proposition [2j The relationship k*(h) in (|5. 10[) is built in accordance with 
(|5.1| . The range of h in <|5 . 1 1 j) is defined by (|5.2|) . The region of admissible values (shaded in Fig. 3) is found using 
the mentioned above fact that for each h —(a + b) the image of Eh under K is a bounded connected segment. 
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// 



Figure 3: The bifurcation diagram of H x K 

Finally, consider the map 

J (4) = H x G : P 6 -> R 2 . 

For < \s\ ^ 6 and \s\ > a, let 
Proposition 5. (i) The bifurcation diagram of the map 

J(4) 

consists of the half-lines Ai, A2, A4, As, A6 and 

the curves 



h = 2s + cj)(s) 
g = — + s 3 + s>(s) 



Ci : <| _ 7 2 , 3 2 , se [-6,0), (5.12) 



s 

h = 2s + <p(s) 

C 2 :{ 7 2 ; * 6(0,6], (5.13) 

# = h« +8 <p(s) 

h = 2s-<j>(s) 

C 3 : { 7 2 3 2 se [a,+oo). (5.14) 

g = — + s -s <f>(8) 

(ii) Denote by go(h) the one-valued function defined by i5.12\) when h ^ —26. Let 

, . _ ( b 2 h-b(a 2 -6 2 ), -(a + 6) < h < -26 
■ 9 * ( j ~\ 9o(/i), Zi^-26 ' (5.15) 

= a 2 /i + a(a 2 -6 2 ). 

TTien i/ie region of admissible values of J*- 4 -* is defined by the inequalities 

g*(h)^g^g*{h), h>-(a + b). 

The bifurcation diagram of J^ 4 ' with the shaded region of admissible values is shown in Fig. 4. 

A straightforward proof of Proposition[5]can be obtained using the same technique as in the proof of Theorem[T] 
Here we just point out some general ideas that explain this result from the point of view of the geometry of S. 

Let v E P 6 be the critical point of J {i ~> . If dH(v) = 0, the n v is an equilibrium, that is, a singular point of the 
system (|3.ip . Such a point is a critical point for each first integral of (|3.ip . In particular, dG(v) — 0. The values 
(h,g) of J^ 4 ) at equilibria (|2.11[) are the points of pairwise intersection of the lines Ai — A4. 

Let 

vaxik{dG(v),dH(v)} = 1, (5.16) 
and c = (g, k, h) = J(v). Then c G S. It follows from (|4.ip - (|4.5[) that (|5 . necessarily implies 

vank{dG(v), dH(v),dK(v)} = 1. (5.17) 
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Figure 4: The bifurcation diagram of H x G 



If the tangent plane to E at the point c is well defined, then the set of zero linear combinations of dG{v), dH(v) 
and dK(v) is one-dimensional. This fact contradicts to (|5 . 1T|) . Therefore, c belongs either to the segment T§ or to 
the set of transversal intersections of two smooth leaves of E. 

The intersection of Ti and r 2 U T3 is nowhere transversal. 

Transversal intersections of Ti and T4 are given by the system (|5.8[) . Solving it with respect to g and h and 
taking into account the admissible region established in Proposition [21 we arrive at the curves (|5 . 1 2|) - (15.14[) . 
Consider intersections of T 2 U T 3 and r 4 . Substitute (|4.4|) for k, g in l|4.7[) : 

(s 2 - a 2 )(s 2 - 6 2 )[2s 2 - 2hs+p 2 } 2 = 0. (5.18) 

Transversal intersections correspond to the values s = ±a, s = ±b (the last multiplier in (|5.18p is responsible 
for the tangency points of r 2 ur 3 and T 4 ). This implies the equations of Ai - A4. As shown in [11) the corresponding 
motions on 01 are the pendulums (|2.12p . (|2 . 13|) . From this fact the inequalities for h are obtained. 

Suppose that T4 has a point of self-intersection. Then for some h the curve defined by (j4.4[) in (g, fc)-plane 
has a double point. Let s\, s 2 be the corresponding values of s. It follows from (|4.4| that 

81 + s 2 = h, s\s\ = 7 2 . 

Hence s%, S2 form one of the pairs ()4.15|) . (|4.16|1 . Substituting these pairs in (|4.4|) gives (I4.14|) . The obtained set of 
points in (h, g) -plane united with the projection of the segment T§ forms the half- lines A5 and Xq. 
The admissible region for (h, g) is established in the same way as in the previous case. 

Propositions [4] and [5] give the explicit formulae (|5.10[) , (|5. 15[) for the values (|5.4p . Then for each h we can 
compute the limits for the parameter s in (|4.4[) corresponding to r 4 n A. Thus the set is completely determined. 
Finally, A/j is obtained as the span of the curves Tj D A^. 

6 Conclusion 

At this point we can draw all bifurcation diagrams of the induced momentum maps on iso-energetic surfaces, 
which are typically five-dimensional. A lot of information on the stability of the critical integral manifolds may be 
immediately obtained for the tori in 971 and 01. The investigation of the new critical set D waits to be fulfilled. 

Since each Eh is a foliation into three-dimensional tori with some degenerations, we can construct the base 
Bh for such a foliation just by factorizing Eh, more exactly, by identifying points of the same connected component 
°f Jg,kh- Then Bh is a two-dimensional analog of Fomenko's graph [S] for an iso-energetic manifold of integrable 
system with two degrees of freedom. In its turn is a bundle over A^ whose fibres are finite sets; the number of 
elements in any fibre is equal to the number of connected components of the corresponding integral manifold. The 



14 



problem of finding this number for all possible situations seems solvable. Then we obtain a complete description 
of the "coverings" Bh — > and, consequently, establish the topology of Bh- 

Naturally, the next step requires new mathematical ideas on how the tori in Eh glue together along the paths 
in the admissible regions. 

If we consider Bh as a two-dimensional cell complex, then, for regular levels of energy, 0-cells correspond 
to closed orbits, a point of each 1-cell represents a two-dimensional torus, and a point of each 2-cell represents 
a three-dimensional torus. The union of the cells of dimensions and 1 forms a graph, to which the method of 
marked molecules [15] can be applied without any modification. The question is what kind of a numeric mark 
should be attached to each two-dimensional cell to obtain from Bh the complete invariant of Liouville's foliation 
of the iso-energetic surface? 

Another approach is to consider the set E° of singular points of (self-intersections, tangency points, and 
cusps), which is easy to obtain from the above results, and associate to each c€ £j the marked loop molecule [T6] . 
In this case, of course, the notion of a mark should be changed to suit increased dimensions of the tori. 

We see that the Kowalevski top in two constant fields provides a highly non-trivial example of integrable 
Hamiltonian system and a complete description of its phase topology is really a challenging problem. 

Received 09.04.05 
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